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Abstract 



o 

CN| In an earlier paper of Cadek, Vokfi'nek, Wagner, and the present au- 

thors, we investigated an algorithmic problem in computational algebraic 
Q topology, namely, the computation of all possible homotopy classes of maps 

^ | between two topological spaces, under suitable restriction on the spaces. 

We aim at showing that, if the dimensions of the considered spaces are 
bounded by a constant, then the computations can be done in polynomial 
time. In this paper we make a significant technical step towards this 
goal: we show that the Eilenberg-MacLane space K(7L, 1), represented as 
U a simplicial group, can be equipped with polynomial-time homology (this 

is a polynomial-time version of effective homology considered in previous 
^ works of the third author and co-workers). 

O To this end, we construct a suitable discrete vector field, in the sense 

of Forman's discrete Morse theory, on K{7L, 1). The construction is purely 
combinatorial and it can be understood as a certain procedure for reducing 
£> finite sequences of integers, without any reference to topology. 

(N The Eilenberg-MacLane spaces are the basic building blocks in a Post- 

nikov system, which is a "layered" representation of a topological space 

\q suitable for homotopy-theoretic computations. Employing the result of 

this paper together with other results on polynomial-time homology, in 
another paper we obtain, for every fixed k, a polynomial-time algorithm 
for computing the fcth homotopy group tt^X) of a given simply connected 

i-H space X, as well as the first k stages of a Postnikov system for X, and 

also a polynomial-time version of the algorithm of Cadek et al. mentioned 

• above. 

X 

1 Introduction 

Recently our co-authors and we |CKM + fl] have developed an algorithm for a 
problem in computational algebraic topology (more precisely, in computational 
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homotopy theory), namely, computing all homotopy classes of maps between 
two topological spaces X and Y (given as finite simplicial complexes, say), 
under certain natural conditions on X and Y. 

Our original motivation was understanding the computational complexity 
of the Zj2-index of a given Zj2-space, which is a quantity appearing in various 
applications of topology in combinatorics and geometry (e.g., topological lower 
bounds for the chromatic number of a graph, or an algorithm for testing the 
embeddability of a given simplicial complex into R rf ). We hope to reach re- 
sults in this direction in the future, and we also expect that the developed 
methods will be applicable for other natural problems (such as extendability of 
maps; as a concrete application, it was already possible to answer a question 
of Franek et al. [FRZ11] on testing nullhomo topy of m aps into a sphere). For 
more information on this project we refer to jCKM+llj . 

Towards polynomial-time homology. The implementation of some of 
the operations in the algorithm of (CKM + 11 relies on the methods of effective 
homology, initiated by the third author in |Ser9 4j and further developed by him 
and his co-workers (see, e.g., [RS02. RRS06, RS12]). These provide algorithmic 
solutions of many problems in algebraic topology, but so far no analysis of their 
running time was available, and for some parts the running time can actually 
be exponential. 

One of our aims is to obtain polynomial-time algorithms for these tasks 
where possible, or alternatively, show computational hardness. 

Let us stress that by "polynomial-time" we mean, throughout this paper, 
polynomial-time for every fixed dimension. Thus, assuming that the input to 
an algorithm is a space represented as a finite simplicial complex X, we want 
that the running time is polynomial in the number of simplices of X, but the 
polynomial may depend on the dimension k of X (and the dependence on k may 
be exponential or even worse). Of course, one could be even more ambitious 
and ask for a polynomial dependence on k as well; however, we do not expect 
such algorithms to exist. 

To integrate this effort with existing algorithms, we start with the framework 
of effective homology mentioned above, and we introduce an analogous defini- 



tion of polynomial-time homology; see Section 2l In another paper CKM + 12 



we show that various known constructions and operations on objects with effec- 
tive homology have polynomial-time versions. With a repertoire of such opera- 
tions, we also obtain a polynomial-time version of the algorithm of |CKM+11 



as well as other algorithms, such as computing the higher homotopy group 
■Kk{X) in polynomial time for every fixed k, or computing the first k stages of 
a Postnikov system for X. 

This paper. Here we make a significant step in this development. First 
we set up the framework of polynomial-time homology (modeled after effective 
homology mentioned above) and some tools of general applicability. Then, 
in the second part of the paper, we present our main technical result. The 
problem which we solve can be formulated purely combinatorially, although in 
this language it perhaps doesn't sound extremely natural: it is a question about 
reducing finite sequences of integers by certain simple operations. We will state 
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it below, and no topological notion at all is required for understanding this 
problem and our solution. 

However, to explain its role in computational topology, we first need to 
sketch some background information. A standard reference for this material 
is May |May67| ; a concise overview is given in |CKM+ll| . and more leisurely 
explanations can be found in [Ser08| or |RS12| . 

A common technique in mathematics and in computer science is to decom- 
pose a general, presumably complicated object into simpler building blocks. For 
the purposes of understanding continuous mappings going into a given topo- 
logical space Y, a suitable decomposition is a Postnikov system for Y; indeed, 



this is a crucial ingredient of the algorithm in CKM + 11 



We do not need to define the rather complicated notion of Postnikov system 
here; it suffices to say that its "building blocks" belong to a particular class of 
topological spaces, called Eilenberg-MacLane spaces and denoted by K(G,k), 
where G is an Abelian group and k > 1 is an integer. In the Eilenberg-MacLane 
spaces appearing in a Postnikov system for Y, the role of the group G is played 
by the homotopy groups 7Tj(Y), i > 2. 

In topology, K(G, k) is defined as a topological space T whose homotopy 
groups satisfy TTk{T) = G and vrj(T) = for alii ^ k. It is determined uniquely 
up to homotopy equivalence (in the class of all CW complexes). 

Generally speaking, the spaces K(G, k) are infinite-dimensional and they 
do not look like very simple objects (with the exception of K(Z, 1), which is 
homotopy equivalent to the circle S 1 ). However, they are in some sense the 
simplest possible spaces concerning maps going into them. These spaces are 
of basic importance in algebraic topology, and a lot of work has been devoted 
to studying their properties, and in particular, computing their homology and 
cohomology (Serre [Ser53] and H. Cartan [Car56j are two of the most famous 
classical works; see, e.g., Clement [Cle02] for an overview and some computa- 
tional aspects). We also refer to Romero and Rubio [RRllJ for an algorithmic 
study of K(G, 1) for noncommutative groups G. 

For the intended algorithmic use, we need a particular representation of 
K(G, k); namely, we need it represented as a particular kind of a simplicial set 
(simplicial sets will be briefly introduced in Section [2] below) , a so-called Kan 
simplicial set. We use the standard Eilenberg-MacLane simplicial model for 
K(G, k); see [EML531 Chapter III], |May67j Chapter V]. 

For the algorithms, we need to equip the simplicial Eilenberg-MacLane 
spaces with polynomial-time homology. The K(G, k) we may encounter can 
have any finitely generated Abelian group as G, and any positive integer as k. 

However, in this paper we will deal only with K(Ti, 1), which serves as a base 
case, while the other K{G, k) can be obtained from it using several operations. 
First, for direct products of groups, we have K(G X H,k) = K(G, k) x K(H, k), 
and so, with a general product operation available, we may assume that G is 
cyclic. Second, a general construction, known as the classifying space (actually, 
in the simplicial setting, we deal with the so-called H^-construction) , allows one 
to pass from K(G,k) to K(G,k + 1), so indeed k = 1 is the important base 
case. Finally, polynomial-time homology for K (Wt/m'Z, 1) can be obtained from 
that for K (Z, 1) using another operation, namely, computing the base space 
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of a fibration. These redu ctions are discussed in [RS12], and polynomial-time 
versions are discussed in |CKM+12j : here we just wanted to provide a quick 



explanation of why the K{7L, 1) case deserves special attention^ 

The combinatorial problem about integer sequences. The /c-dimensional 
simplices of the standard simplicial model of KiTL, 1), k = 0, 1, . . . can be rep- 
resented by /c-term sequences of integers. With the traditional "bar notation" , 
such a sequence is written as 

o- = [oi | a 2 | • • • | Ofc], ai,a 2 ,...,a k eZ. (1) 

In the rest of this introduction, a "^-dimensional simplex" will thus be synony- 
mous with a "A:-term sequence of integers" . 

For our problem we consider only nondegenerate simplices, represented by 
sequences with no zero terms. Thus, from now on, we always assume that all 
the dj are nonzero. 

For each k, there are k + 1 face operators do, d%, . . . , d k , which map /c-term 
sequences to k — 1 term sequences: do deletes the first component, d k deletes the 
last component, and for i = 1,2, . . , ,k— 1, <9j reduces the number of components 
by one by adding together the ith. and (i + l)st component. More formally, with 
a as above, 

d a = [a 2 | • • • | a k ], d k a = [ai \ ■ ■ ■ \ a k -i], 

did = [a\ | • • • | di-i | ai + a i+ i | a i+2 \ ■ ■ ■ \ a k ], 1 < i < k — 1. 

The goal is to divide the set of all possible finite sequences a of nonzero 
integers into three classes S, T, and C (the source simplices, target simplices, 
and critical simplices), and construct a bijection V: S — > T (which will be 
called a discrete vector field), such that for every a £ S, we have a = d{V(a) 
for exactly one i. We also require certain additional properties, which we explain 
next. 

With S,T,C, and V as above, let us consider a sequence (simplex) a £ S 
of some dimension k, and let us say that a simplex r (of dimension k or k + 1) 
is reachable from a if it can be reached from a by finitely many moves, where 
the allowed moves are 

• passing from a current simplex a G S to the simplex r = V(a) G T, and 

• passing from a current simplex r G T to a simplex a = c^t G 5 U C such 
that r 7^ V{a), where i G {0, 1, . . . , k). 

With these definitions, it is required that 



1 Curiously, K{7L, 1) as a topological space almost can't be simpler — as we mentioned, 
it is homotopy equivalent to the circle S , and other Eilenberg-MacLane spaces are much 
more complicated. But we need to work with the Kan simplicial model of K{7L, 1) as intro- 
duced above, which has infinitely many simplices in every dimension k > 1. As we will see, 
for effective (or polynomial-time) homology, it is not sufficient to know, for example, that 
H2(K(7L, 1)) = 0, but we need to be able to actually compute "witnesses" for it; that is, given 
a 2-cycle Z2 on K(7L, 1), compute a 3-chain for which 22 is its boundary. This problem would 
be trivial for the standard simplicial representation of S 1 with one vertex and one edge, but 
it is not trivial for the considered Kan model of K(7L, 1). 
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(i) for every k, C contains only finitely many fc-dimensional simplices; and 



(ii) starting with any a, we can never make an infinite sequence of allowed 
moves; that is, we can reach only finitely many simplices, and we also 
cannot get into a cycle. 

Moreover, we measure the size of a simplex a = [a± \ ■ ■ ■ |a^] as the total num- 
ber of bits needed to write down ai, . . . , a^; more formally, we set size(cr) := 
Yli=i size(aj) and size(a) := 1 + Llog 2 (H + 1)J . Then we also require that 

(iii) For every /c-dimensional simplex a, the sum of size(c) over all a reachable 
from a is bounded by a polynomial (depending on k) in size(a). 

To illustrate these definitions, let us present a classical vector field Veml 
due to Eilenberg and Mac Lane, which satisfies (i) and (ii) (and yields effective 
homology for K{Ti, 1)) but not (iii). 

There are only two critical simplices, the O-dimensional [] (the empty se- 
quence) and the 1-dimensional [1] {^] The set S of source simplices consists of 
the sequences with a\ ^ 1, while T contains the sequences with a± = 1 (the 
two critical simplices are exceptions to this rule). 

For a = [ai\ • • • ja^] S 5, ai ^ 1, the vector field Veml is defined by 



It can be checked that, for any starting a, the sequence of moves is determined 
uniquely (there is no branching). We will not verify this, but it may be instruc- 
tive for the reader to write down small examples of sequences and see how the 
moves work. 

It is easy to see that, for a positive integer a, the sequence of moves starting 
from [a] is [a] — > [l\a — 1] — > [a — 1] — > [l\a — 2] — > [a — 2] — > . . .; there are about 
a moves, and this is exponential in the number of bits of a. Thus, condition 
(iii) above indeed fails. 

We will provide a solution satisfying (i)-(iii) in Section |4j Before that, we 
introduce simplicial sets, polynomial-time homology, and discrete vector fields 
in general. 

2 Simplicial sets with polynomial-time homology 

Simplicial sets. A simplicial complex is a way of specifying a topological 
space in purely combinatorial terms, and also a way of presenting a topological 
space as an input to an algorithm; we assume that the reader is basically familiar 
with this concept. 

A simplicial set can be regarded as a generalization of a simplicial complex; 
it is more complicated, but more powerful and flexible. The algorithms we 

2 This actually corresponds to the topological fact that the considered K(7L, 1), as a topo- 
logical space, is homotopy equivalent to S ; [] represents a vertex, and [1] an edge glued to 
that vertex by both ends, forming an S . 




[l|ai — 1 1 a.2 1 • • • for a\ > 1, 
[1 1 ai | a.2 1 • • • \ cik] for a\ < 0. 
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consider use simplicial sets as the main data type for representing topological 
spaces and their maps. A friendly introduction to simplicial sets is [Frillj . and 
another introductory treatment can be found in [Ser08| ; older compact sources 
are, e.g., |Cur7H May67| , and [GJ99j is a more modern and comprehensive 
treatment. 

Similar to a simplicial complex, a simplicial set is a space built of vertices, 
edges, triangles, and higher-dimensional simplices, but simplices are allowed to 
be glued to each other and to themselves in more general ways. For example, one 
may have several 1-dimensional simplices connecting the same pair of vertices, 
a 1-simplex forming a loop, two edges of a 2-simplex identified to create a cone, 
or the boundary of a 2-simplex all contracted to a single vertex, forming an S 2 . 




Another new feature of a simplicial set, in comparison with a simplicial 
complex, is the presence of degenerate simplices. For example, the edges of 
the triangle with a contracted boundary (in the last example above) do not 
disappear, but each of them becomes a degenerate 1-simplex. 

A simplicial set X is represented as a sequence (Xq, X%, X2, ■ ■ .) of mutually 
disjoint sets, where the elements of X^ are called the k-simplices of X (we note 
that, unlike for simplicial complexes, a simplex in a simplicial set need not be 
determined by the set of its vertices; indeed, there can be many simplices with 
the same vertex set). For every k > 1, there are k+1 mappings do, ■ ■ ■ , df~ : X^ — > 
Xk-i called face operators] the intuitive meaning is that for a simplex a G X^, 
diO is the face of a opposite to the ith vertex. Moreover, there are k + 1 
mappings so,...,Sk'- Xk — > Xk+i (opposite direction) called the degeneracy 
operators] the approximate meaning of SjU is the degenerate simplex which is 
geometrically identical to a, but with the ith vertex duplicated. A simplex is 
called degenerate if it lies in the image of some Si] otherwise, it is nondegenerate. 
We write A ndg for the set of all nondegenerate simplices of X. 

There are natural axioms that the <9j and the Si have to satisfy, but we 
will not list them here, since we won't really use them. Moreover, the usual 
definition of simplicial sets uses the language of category theory and is very 
elegant and concise; see, e.g., [GJ991 Section 1.1]. 

Every simplicial set X specifies a topological space \X\, the geometric real- 
ization of X. It is obtained by assigning a geometric /c-dimensional simplex to 
each nondegenerate fc-simplex of X, and then gluing these simplices together 
according to the face operators; we refer to the literature for the precise defini- 
tion. 

There is a canonical way of converting a simplicial complex to a simplicial 
set; basically, one just needs to add appropriate degenerate simplices. 

We have already given a relatively sophisticated example of a simplicial set, 
namely, K(7L, 1), or more precisely, the standard Eilenberg-MacLane represen- 
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tation of K(Z, 1) as a Kan simplicial selj^jas defined in the introduction (except 
that we haven't yet specified the degeneracy operators, which are very simple: 
Si inserts after the ith component of a sequence). 

Representing infinite simplicial sets. In many areas where computer 
scientists seek efficient algorithms, both the input objects and intermediate 
results in the algorithms are finite, and they can be explicitly represented in 
the computer memory; this is the case, e.g., for algorithms dealing with graphs 
or with matrices. 

In contrast, in the algorithms for homotopy-theoretic questions considered 
here and in related works, we need to deal with infinite objects. For example, 
even if the input is a finite simplicial complex, its Postnikov system (mentioned 
in the introduction) is made of Eilenberg-MacLane spaces, such as K(Z, 1), 
represented as Kan simplicial sets, and these are necessarily infinite. More 
concretely, as we have seen, K(Z, 1) has infinitely many simplices in each di- 
mension k > 1 (and thus we cannot explicitly store even the part up to some 
fixed dimension). 

For algorithmic purposes, we thus represent a simplicial set X by a collection 
of several algorithms, which allow us to access certain information about X, 
without having all of it explicitly stored in memory. (In computer science, this 
is also called a black box or oracle representation of X, and in the terminology 
of object-oriented programming, we can think of X as an instance of a class 
"simplicial set".) A similar representation is used for other kinds of infinite 
topological or algebraic objects as well. 

Locally effective simplicial sets. For some computations, it may be suffi- 
cient to represent X by a black box providing only "local" information about 
X, and in that case, in accordance with the terminology in earlier papers, e.g., 
[RS02, RS12. iRSllj . we speak of a locally effective representation. 

Concretely, let X be a simplicial set, and suppose that some computer rep- 
resentation ("encoding") for the simplices of X has been fixed (e.g., in the case 
of K(% : 1), we can fix the representation of the simplices of K (Z, 1) by integer 
sequences, and represent the integers in the sequences by the standard binary 
encoding). We say that X is a locally effective simplicial set if algorithms are 
available that, given (an encoding of) a A;-simplex a of X and i 6 {0, 1, . . . , k}, 
computes the simplex <9j<r, and similarly for the degeneracy operators Sj. Briefly 
speaking, the face and degeneracy operators should be computable maps. 

Computing global information. Suppose that we want to compute some 
"global" information about a given simplicial set X, for example, the kth ho- 
mology group Hk{X). Then a locally effective representation of X is typically 
insufficient, and we need to augment it in some way. 

Of course, in the particular example with the homology groups, we could 
insist that X be augmented with a black box that, given k, returns some rep- 

3 We won't define a Kan simplicial set, but we just mention a key property, which is the 
reason why these simplicial sets are essential to the considered algorithms. Namely, if X 
is a simplicial set and Y is a Kan simplicial set, then every continuous map |X| — > \Y\ is 
homotopic to a simplicial map X — > Y . Thus, continuous maps into Y have a combinatorial 
representation, describing them up to homotopy. 
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resentation of Hk(X). The problem is that X may not be given to us directly; 
rather, we may need to construct it from other simplicial sets by a sequence 
of various operations. For example, in the introduction we mentioned that the 
Eilenberg-MacLane spaces K{G,k) can be constructed starting with K(Ti, 1) 
and applying operations of several kinds, such as product or classifying spaceQ 
Then, for example, a black box for computing the homology groups of X is 
not in itself sufficient to compute the homology groups of the classifying space 



The third author and his co-authors have developed a more sophisticated 
way of augmenting a locally effective simplicial set X with homological informa- 
tion, which is captured in the notion of a simplicial set with effective homology. 
These simplicial sets do possess a black box for computing homology groups, but 
they are also equipped with additional information, which makes them stable 
under a large repertoire of operations: if we apply some of the "classical" opera- 
tions, such as product, classifying space, loop space, etc. to simplicial sets with 
effective homology, the result is again a simplicial set with effective homology 
(and in particular, it has a black box for computing homology groups) F] 




It may be useful to keep in mind that, since a simplicial set is represented 
by a black box, operations on such simplicial sets are performed by composition 
of algorithms; i.e., the black box for the new simplicial set operates by calling 
the black boxes of the old sets and processing the values returned by themj^] 

For defining a simplicial set with effective homology, and their polynomial- 
time counterpart, we need to recall some notions concerning chain complexes. 

Chain complexes. For our purposes, a chain complex C* is a sequence 
(Cfc)fcL_oo of free Z-modules (i.e., free Abelian groups), together with a sequence 
(dk : Ck — » Ck-i)'£L_ 00 of group homomorphismsj^] The Ck are the chain groups, 
their elements are called k-chains, and the dk the differentials. The differentials 
have to satisfy dk-idk = for every h (here dk-\dk denotes the composition 
of maps). We also recall that the fcth homology group Hk(C*) of the chain 
complex C* is defined as the factor-group ker dk/ im4+i. 

For every simplicial set X, there is a canonically associated chain complex, 
which is used to define the homology groups Hk(X). Actually, there are two 
natural possibilities, depending on whether degenerate simplices are taken into 
account. We use the normalized chain complex, which is based solely on the 

As another, perhaps more sophisticated example, we can mention the computation of 
the homotopy group nk(X) for a 1-connected simplicial set X: for this, given X, one first 
produces another simplicial set X' from X, by a sequence of operations that "kill" the first k— 1 
homotopy groups, and then -Kk(X) is computed as Hk(X') using the Hurewicz isomorphism. 

5 One can also consider other kinds of objects with effective homology, such as chain com- 
plexes, but for concreteness, we will stick to simplicial sets. 

6 This feature makes it very natural to implement algorithms from this area using functional 
programming languages, as was done for the package Kenzo; see, e.g., [HPRSlT] . 

7 These chain complexes are over 7L; more generally, one considers chain complexes over 
a commutative ring R, where the Ck are _R-modules. These are needed, among others, for 
homology with coefficients in R. But for our purposes, homology with integer coefficients 
suffices; if needed, homology groups with other coefficients can be computed using universal 
coefficient theorems. Alternatively, all of the theory can be built with coefficients from a fixed 
ring 7?, provided that R is equipped with sufficiently strong algorithmic primitives. 



of X. 
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nondegenerate simplices. We reserve the simple notation C*(X) for it. 

Thus, Cfc(X) denotes the free Abelian group over g , the set of all k- 
dimensional nondegenerate simplices (in particular, Ck{X) = for k < 0). 
This means that a A:-chain is a formal sum 

c = ^2 "cr • cr, 

where the a a are integers, only finitely many of them nonzero. The differentials 
are defined in a standard way using the face operators: for /c-chains of the form 
1 • a, which constitute a basis of Ck(X), we set c4(l • a) := Yli=o(~ ' ®i a 
(some of the diO may be degenerate simplices; then they are ignored in the 
sum), and this extends to a homomorphism in a unique way ("linearly"). 

We note that if X is a locally effective simplicial set, then the fc-chains 
of C*(X) are finite objects; a /c-chain c can be represented by a list of the 
fc-simplices a on which c is nonzero, and of the corresponding coefficients a a . 
Then the differentials are computable maps. 

However, if X^ ds is infinite, then Ck(X) has infinite rank, and we cannot 
use it directly for computing homology groups. The solution adopted in ef- 
fective homology is to have, together with a locally effective simplicial set X, 
a reduction from C*(X) to an "effective" chain complex EC*, for which each 
chain group EC^ has a finite rank. 

Reductions. Let C*, C* be two chain complexes. To define a reduction from 
to , we first recall two other standard notions from homological algebra: A 
chain map f : C* — > C* is a sequence (/fc)feL-oo of homomorphisms f^: Cy. — )■ 
compatible with the differentials, i.e., fk-idk = dkfk- If f,9' C* — > (7* are two 
chain maps, then a chain homotopy of / and g is a sequence (fyfc)fcL— re °^ 
homomorphisms h^: — > Cfc+i such that f — g = d^+ih^ + h^-idk- 

Now a reduction p from C* to C* consists of three maps /, g, h, such that 

• / : C* — > C* and <? : C* — >• C* are chain maps; 

• the composition fg: (7* — > (7* is equal to the identity id^ , while the 
composition gf : C* — >• C* is chain-homotopic to idc^ , with h : C* — > C* 
providing the chain homotopy; and 

• fh = 0,hg = 0, and /i/i = 0. 

The notion of reduction goes back to Eilenberg and Mac Lane [EML53| 
Section 12], who called it a contraction^] It is routine to check that if there is 
a reduction from C* to (7*, then C* and C* have isomorphic homology groups 
in each dimension. Reductions can also be composed, as follows: if (f,g,h) is 

a reduction from C* to (7* and (f',g',h') is a reduction from C* to (7*, then 

(/'/, gg' , h + gh! /) is a reduction from (7* to C*. 

8 They did not require the condition hh — 0, but simple transformation converts a reduction 
without this condition into another one satisfying it. 
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Effective homology. We are getting close to stating the definition of a 
simplicial set with effective homology. The last step is to define what we mean 
by an effective chain complex EC*. We assume that, first, EC* is locally effec- 
tive, meaning that each chain group ECk has some distinguished basis Bas^, 
&:-chains are represented as linear combinations of elements of Basfc (and thus 
they can be added and subtracted algorithmically) , and there is an algorithm 
for evaluating the differentials df.- Second, EC* is effective, which means, in 
addition to the above, that there is an algorithm that, given k, outputs the list 
of elements of the distinguished basis Bas^; in particular, this implies that each 
ECk has a finite rank ry.. We note that by combining the construction of Basfc 
with the ability to evaluate the differential dk, we can compute the matrix of 
dk with respect to the distinguished bases Bas^ and Basfc_i. 

We can now define a simplicial set with effective homology as a locally ef- 
fective simplicial set X together with an effective chain complex EC* and a 
reduction p from C*{X) to EC*, where the three maps f,g,h from the defini- 
tion of reduction are computable^] 

In this paper we won't have the opportunity to demonstrate the usefulness of 
effective homology in algorithms; we refer to, e.g., |RS12l SerOB] for examples of 
applications. The following short remark may provide at least some illustration 
of how things work. 

Let us assume that we are (explicitly) given a /c-cycle c G Ck(X), where 
X is a simplicial set with effective homology as above, and we would like to 
test whether it is its boundary, i.e., whether c = dk+ib for some b G Ck+i(X). 
We compute c := fk(c) G ECk, and test if it is a boundary in EC*. If not, 
then neither c is a boundary. Otherwise, we can compute b G ECk+i with 
c = dk + \b. Then, as is easy to check, c is the boundary of b := g(b) + h(c). 
In this sense, the reduction provides not only a way of recognizing boundaries, 
but also computable witnesses of being a boundary. 

Polynomial-time homology. The meaning of polynomial-time homology for 
the simplicial set K{%, 1) considered in this paper is defined in a straightforward 
way: we want the face and degeneracy operators to be computable in polynomial 
time (which is obvious in this particular case), and K (Z, 1) should be equipped 
with effective homology as above in such a way that, for every k, the maps 
fk>9k,hk are computable in polynomial time, with the polynomial possibly 
depending on k as usual. 

We stress that since we deal with a single effective chain complex EC*, the 
ranks r& depend only on k and thus, for k fixed, they are constants. The matrix 
of the differential dk in EC*, too, is a constant-size object. 

However, our setting with K(Z, 1) is somewhat unusual in the analysis of 
algorithms: We are dealing with a single simplicial set, fixed once and for all, 
which does not depend on any input. This is an exceptional setting; most algo- 
rithms work with objects that do depend on the input. To draw an analogy from 

9 In [RSI 2] and in other papers, effective homoiogy is defined in a more generai way, using 
strong equivalence of chain compiexes instead of just a reduction. A strong equivalence of 
C* and C* means that there is an auxiliary chain complex A* and reductions of A* to both 
C» and C*. However, here the simpler notion using a single reduction suffices, and this only 
makes the result formally stronger, since a reduction is a special case of a strong equivalence. 
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a different area, the setting of the present paper can be compared to seeking an 
algorithm for computing the nth digit of the number tt, while the more usual 
case would be to consider algorithms for evaluating arithmetic expressions with 
arbitrary precision, where we start with integer numbers as inputs and apply 
addition, subtraction, multiplication, division, roots and functions like exp, In 
or arcsin. 

To have an example from the area considered here, in an algorithm for 
computing with a given topological space X, say specified as a finite simplicial 
complex, we may need polynomial-time homology for the Eilenberg-MacLane 
space K{W l , 1), where n is a parameter depending on X. Then we want that 
in the corresponding effective chain complex for K(Ti n , 1), the ranks r2,r%, etc. 
each depend polynomially on n. (Of course, for this to be useful, we also need 
that n depends at most polynomially on the size of X.) 

This example suggests that, in order to have a generally useful notion of 
polynomial-time homology, we need to define it formally for a whole family, 
ty pically infin ite, of simplicial sets. Here we present this issue briefly, referring 
to CKM + 12 for a more detailed discussion. 

Let X be a set, typically countable, such that each element I £ I has some 
agreed-upon computer representation (i.e. encoding by a finite string of bits). 
A simplicial set parameterized by I is a mapping X that assigns a simplicial 
set X{I) to each I £ I. We also assume that the simplices of each X{I) 
have some encoding by bit strings. Then we define a locally polynomial-time 
simplicial set as a simplicial set X parameterized by some X such that the face 
and degeneracy operators on a fc-simplex a of X(T) can be evaluated in time 
polynomial in size(J) + size(cr), where the polynomial may depend on k (and 
size(.) denotes the the number of bits in the encoding). 

Quite analogously, we define a chain complex C* = (C(J)* : I € I) param- 
eterized by a set X. We say that such a C* is locally polynomial-time if each 
C(i% is a locally effective chain complex (and in particular, it has a distin- 
guished basis Bas(J)fe, and /c-chains are represented w.r.t. this basis), and for 
each fixed k, the differential (c£j)& on C(I)f. can be evaluated in time polyno- 
mial in size(J) plus the size of the input fc-chain. We note that addition and 
subtraction of /c-chains are polynomial-time operations automatically. 

We say that a simplicial set X parameterized by a set X is equipped with 
polynomial-time homology if the following hold. 

• X is locally polynomial-time. 

• There is a locally polynomial-time chain complex EC* , also parameterized 
by X, such that, for each fixed k, the distinguished basis Bas(/)fc of EC{I)k 
can be computed in time polynomial in size(J), and in particular, the rank 
r(I)fc is bounded by such a polynomial. 

• For every I £ I, there is a reduction pi from C*{X{I)) to EC(I)*, where 
the maps (//)&, (<7/)fc, (hj)k of pi are all computable in time bounded by 
a polynomial in size(I) plus the size of the input /c-chain; the polynomial 
may depend on k. 
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3 Polynomial- time homology from a discrete vector 
field 



Discrete Morse theory, developed by Forman [For98 (also see |For02j), belongs 
among fundamental tools in combinatorial topology. For us, the important 



point is that a suitable discrete vector field on a simplicial set 10 X can be used 
to equip X with effective homology; this is an implication of one of Forman's 
results, as was observed by Romero and Sergeraert |RSllj (they also generalized 
Forman's construction by dropping a certain finiteness condition). Here we 
review the definitions, more or less repeating in a general setting the definitions 
given for K(%, 1) in the introduction. Then we formulate a sufficient condition 
on the vector field so that the construction provides polynomial-time homology 
for X. 

Discrete vector fields. Let X be a simplicial set. For a simplex r G X, 
it may happen that two face operators give the same simplex, i.e., diT = djT, 
i 7^ 3 (geometrically, this means that the two faces of the simplex r are "glued 
together"). We say that a is a regular face of r if a = d{T for exactly one 
index i. 

A discrete vector field V on a simplicial set X is a set of ordered pairs 
(directed edges) of the form (a, r), where a, r G X ndg , a is a regular face of r, 
and for every two distinct pairs (a, r), (a' , t') G V, all of a, r, a' , r' are distinct. 

Given a discrete vector field V, the nondegenerate simplices of X are clas- 
sified into three subsets S, T, and C as follows: 

• S are the source simplices; these are simplices a such that (<r, r) G V for 
some t. 

• T are the target simplices; these are simplices r such that (a, r) G V for 
some a. 

• C are the critical simplices; these are the remaining simplices, not occur- 
ring in any edge of V. 

Often it is useful to regard V as a bijective mapping V: S — > T, as we did 
in the introduction. Thus, for (a, r) G V, we sometimes write r = V(o~) and 
<j = V- 1 {t). 

In a drawing of a simplicial set, the pairs (a, r) of a vector field can be 
indicated by arrows pointing from a into r, as in Fig. [TJ 

Admissible vector fields and the V9-graph. The vector fields useful in 
discrete Morse theory, as well as in our context, have an extra property. For 
defining it, we first introduce an auxiliary directed graph, as drawn in Fig. [2j 
which we call the V9-graph. 

The vertex set of the Vd- graph is X ndg . In the drawing, the empty cir- 
cles correspond to source simplices, the full circles to target simplices, and the 
critical simplices are marked by double circles. 



10 In |RS11| . vector fields are considered in somewhat greater generality, on algebraic cell 
complexes. Here it is sufficient to stay in the perhaps more intuitive setting of vector fields on 
simplicial sets. 
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Figure 1: A triangulation of the real projective plane with a discrete vector field 
(after Forman [For02j, Fig. 4.1). Pairs of vertices with the same label should 
be identified; thus, there are only one critical edge and one critical vertex. 




The edges of the Vd- graph are of two kinds: first, those belonging to V 
(drawn bold and pointing upwards), and second, all edges of the form (r, a), 
where r is a target simplex, a is a face of r and a source or critical simplex, 
and (cr, r) V (these edges point downwards) j^j These edges correspond to 
the "allowed moves" defined in the introduction. 

We call the vector field V admissible if the V9-graph contains no directed 
cycle and no infinite directed path. The field in Fig.[T]is admissible, for example. 

One of Forman's results says that an admissible vector field V can be used to 
"simplify" the underlying simplicial set X: by a sequence of suitable collapsing 
operations, which is defined based on V, one obtains a cell complex (no longer 
necessarily a simplicial set), which is homotopy equivalent to X but typically 
much smaller — its cells correspond only to the critical simplices. 

We will not use this result directly (and thus we don't formulate it precisely). 
Rather, we build on a related result (obtained implicitly by Forman with an 
additional finiteness assumption, and explicitly and in general in |RSllj ). as- 
serting that an admissible vector field provides a reduction of the normalized 
chain complex C*{X) to a suitable chain complex C£ nt . In this chain complex, 
each C£ nt is the free Abelian group on the set of all A;- dimensional critical sim- 
plices. The differentials in C£ nt are defined based on V, and they are locally 
effective assuming that X and V are locally effective in a natural sense. 

Polynomially bounded vector fields. We need a polynomial-time version 
of this result. Let V be an admissible vector field V on a locally polynomial- 
time simplicial set; we assume that both X and V are parameterized by a set X, 
as in the definition of a locally polynomial-time simplicial set 
let reachy (c) (or just reach (cr) if V is understood) denote the set of all simplices 
reachable from cr by a directed path in the V9-graph. 

Let us say that V is polynomially bounded if the following hold: 

(PBV1) An algorithm is available that, given I £ I and a simplex a G X(I)^ d& , 
classifies a as source, target, or critical. In the source case, it also returns 
the simplex Via). The running time is polynomial in size(/) + size(c) for 
every fixed k. 

(PBV2) For every fixed k and every a G X(I)^ dg , the sum of encoding sizes of all 
simplices in reach v(°~) is bounded by a polynomial in size(J) + size(cr). 

Theorem 3.1. If X is a locally polynomial-time simplicial set and V is a 
polynomially bounded vector field on X such that, for every k, the sum of the 
encoding sizes of all k-dimensional critical simplices is polynomially bounded 
(in s\ze(I)), then X can be turned into a simplicial set with polynomial-time 
homology. 

In a simplicial set, it may happen that cr is a "multiple" face of r. i.e., a = diT holds for 
several indices i. In such case, we connect r to a with multiple edges in the V9-graph, one 
edge for each such index i. 

12 Of course, for the main result of this paper, polynomial-time homology for K(%, 1), pa- 
rameterization is not needed, but we need it if we want to have a general tool for obtaining 
polynomial-time homology from a vector field. 



For a G X nd §, 
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Proof. The proof essentially follows by inspecting the work of Forman [For98 
(mainly Sections 7 and 8) and making simple observations about the compu- 
tation of the relevant maps. For the reader's convenience, we provide a self- 
contained presentation (this seems simpler and not much longer than referring 
to the appropriate claims in Forman's paper, introducing his notation, etc.). 
Our presentation is, similar to that of Forman, mainly in a combinatorial lan- 
guage. We refer to [RSllJ for two other, more algebraic variants of essentially 
the same proof. 

Throughout the proof, we keep the parameterization of X and V by I 
implicit. 

To provide the desired reduction from C* := C*(X), we need to define the 
target chain complex C£ nt and provide the three maps /, g, h as in the definition 
of a reduction. We begin with introducing several auxiliary maps and checking 
some of their properties. 

The vector field V induces a sequence V# = (V#k)kL-oc °^ homomorphisms 
V#k- Ck Cfc+i, as follows: for a source fc-simplex a, we have V#fc(l • a) := 
(— • V(a), where i is the unique index with a = diV(a), and for a target 
or critical, we have V#fc(l • a) := 0. 

Next, we introduce a chain map f> : C* — > C* by 

$ := 1 + V#d + dV#, 

where 1 stands for the identity chain map and d is the differential of C*. It 
is easy to check that $ is a chain map: indeed, d& = d + dV#d + ddV# = 
d + dV#d = $d (using dd = 0). 

For the proof, it is important to understand how <3> works. We will thus 
discuss how the image 3>(l-cr) is formed, depending on the type of a fc-simplex a. 

1. The simplest case is a a target simplex; see Fig.[3]left. Then V#(l-a) = 0, 
and thus $(1 ■ a) = 1 • a + Yli=o 1)* ' di a )- So we consider all 
faces a' of a, with the appropriate signs, and apply V# to them. Only 
the a 1 that are sources may contribute to the image (and then (a, a') are 
edges of the V9-graph), and <3?(1 - a) is supported only on target simplices. 

Moreover, we observe that, crucially, the coefficient of a in $(1 • a) is 
0; indeed, if j is the unique index with V~ l (a) = djO~, then we have 
V # ({-l) j ■ dja) = (-iy +1 (-iy ■ a = -1 • a, which cancels out with 
the 1 • a coming from the 1 in the definition of (Here we rely on 
the condition that V~ l (o~) is a regular face of a from the definition of 
discrete vector field, since we need the coefficient of V^ 1 (a) in d(l • a) to 
be invertible, i.e., equal to ±1.) 

Summarizing, $(1 • a) consists of the target simplices reachable from a in 
exactly two steps in the V9-graph, with appropriate signs. 

2. For a a critical simplex we find, by a similar reasoning, that $(1 • a) 
consists of a with coefficient 1, plus all the (target) simplices reachable 
from a in exactly two steps in the V9-graph, again with appropriate signs. 
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Figure 3: Forming the image $(1 • a). 

3. Finally, for a a source, both the dV# and V#d terms may make a nonzero 
contribution to <3?(1 • a). For dV# (going first up, then down), we get, with 
appropriate signs, all the source simplices reachable from a in exactly two 
steps in the V9-graph, with a itself cancelled out, plus some additional 
target and critical simplices (here we do not follow the edges of the Vd- 
graph — that's why the arrows are dotted in the picture). For V#d (first 
down, then up), we get only target simplices. 

Next, we define = limjv-»oo $> N as the stabilization of <£; that is, given 
a A;-chain c, we compute <3>(c), ^(^(c)), etc., until we reach a chain c with 
3>(c) = c, and we set <£°°(c) := c. 

To check that the iterations of indeed stabilize after finitely many steps, 
it suffices to consider the case c = 1 • a, and then the stabilization follows easily 
from the above discussion of the action of (and from the admissibility of the 
vector field V). Moreover, we can see that the chains in im$°° are supported 
only on critical and target simplices. 

We also need to check that <fr°° is computable in polynomial time. In or- 
der to compute <£°°(1 • a) (which is sufficient), we just compute the iterations 
$ N (1 ■ a), N = 1, 2, . . ., until they stabilize. We note that each simplex in the 
support of some $> N (1 ■ a) can be reached from a by following a directed path 
in the V9-graph, then possibly going to a face of the current simplex (a step 
corresponding to a dotted arrow in Fig. [3]), and then again following a directed 
path in the V9-graph. Hence, by the polynomial boundedness of the vector 
field V, the stabilization occurs for iV at most polynomially large, and the sum 
of the encoding sizes of all simplices in the supports of all chains encountered 
along the way is also polynomially bounded (essentially by the square of the 
bound in condition (PBV2)). 

Each coefficient in the chain $ Ar+1 (l • a) is the sum of O(k) coefficients in 
$> N (1 - a). So each coefficient in $> N (1 ■ a) is bounded by exp(0(iV)), and hence 
its size (number of bits) is at most O(N). Therefore, $°° is indeed polynomial- 
time computable. 



supp $(1 • a) 
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Now we define an auxiliary chain complex Cf; we set C* := im$^° C Ck- 
Equivalently, as is easily seen, C* = {c 6 Ct : <fr(c) = c}. The differential of C* 
is the restriction of the differential of C* (this works since $ is a chain map). 
Let i: C* — > C* be the inclusion (which is a chain map). 

Next, we come to the definition of C£ nt ; as was announced above, the 
chain group Cf is the free Abelian group (Z-module) with the set of the 
A;-dimensional critical simplices in X as a basis. It remains to define the differ- 
ential. 

First we let jk '■ C* — > C£ rit be the homomorphism that restricts a chain 
c G C* to the critical simplices (i.e., for c = ^2 ae x k a a ■ cr, we set jk(c) := 
J2aex k r\C a a' a )- We observe that viewed as a homomorphism C£ nt — > C*, 
is an inverse to jk- Indeed, from the description of <& given above, it is easy to 
see that for a critical, <j?°°(l - a) = 1 • a + d for some c' supported on target 
simplices, and from this the claim follows. 

Hence each C£ nt is isomorphic to C* , and the differential d CTlt of C£ rit can be 
defined so as to make j and mutually inverse chain isomorphisms; explicitly, 
^cnt ._ This finishes the definition of the target chain complex for 

the desired reduction; it is clear that the matrices of the differential d cnt are 
polynomially computable, provided that the total encoding size of the critical 
simplices is polynomial in each dimension. 

It remains to define the maps /, g, h in the reduction. The following diagram 
summarizes the relevant chain complexes and maps defined so far, plus /, g, h: 



f 




fJ 



As the diagram suggests, we put / := j<3>°° and g := Then, since j and 

<£°° are mutually inverse and = 1, we have fg = 1, as required by the 
definition of a reduction, and gf = i<&°°. 

The chain homotopy h of i<&°° with the identity (Forman uses the letter L 
for this map) is now defined as the stabilization of the maps 

-y # (l + $ + $ 2 + --- + $ N ), N = l,2,... 

To see that these iterations indeed stabilize on each chain 1 -a, we recall that for 
sufficiently large N, & N (1 ■ a) is supported only on critical and target simplices, 
and V# sends such chains to 0. By essentially the same argument as that for 
the computability of we also get that each hk is computable in polynomial 
time. 

We need to verify that h is the required chain homotopy, i.e., dh + hd = 
1 — This is a simple formal calculation (showing where the formula for h 

comes from), which we leave to the reader (also see [For98, proof of Th. 7.3]). 

As the last step, we want to check the conditions fh = 0, hg = 0, and 
hh = 0. To this end, we note that the chains in im/i are supported only on 
target simplices. Moreover, if c is a chain supported only on target and critical 
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simplices, then $(c) has the same property, and hence h(c) = 0. These two 
properties immediately give hh = 0. Similarly, img = im$°° is supported only 
on target and critical simplices, and hence hg = 0. Finally, we have seen that 
maps target simplices to 0, and so does / = , which gives fh = and 



concludes the proof of Theorem 3.1 □ 



4 A polynomially bounded vector field for K{1L 1 1) 

Here we finally get to the combinatorial core of the paper; we will provide a 
polynomially bounded vector field for KiTL, 1). 

A simple composition of vector fields. For the sake of presentation, it 
will be easier to split the vector field into two parts. Roughly speaking, the 
first part will get rid of all negative components in the considered sequences 
\a\ \ ■ ■ ■ |ofc], and the second part will do the rest. 

Here is the way of "splitting into two parts" in a general setting. Let X be a 
simplicial set, let V\ be a vector field on X, with the set C\ of critical simplices, 
and suppose that C\ is closed under the face operators (each face of a critical 
simplex is again critical, or degenerate). Let Y be the simplicial subset of X 
induced by C\ (i.e., its nondegenerate simplices are the critical simplices of Vi), 
and let V2 be a vector field on Y. 

Then we can define a "composition" V of V\ and V2 in the obvious way; 
formally, if we regard a vector field a set of ordered pairs, we simply set V := 
V\ U V2. Clearly, V is a vector field, and it is easily seen that V±, V2 admissible 
imply V admissible, and similarly for polynomial boundedness. 

In the case of X = K(Z, 1), the role of Y will be played by the simplicial 
set whose simplices are the integer sequences with all terms nonnegative. With 
some abuse of the usual notation, we will denote this simplicial set by K(N, 1). 

The first vector field will be denote by Vb s and called the bubblesort field, 
since directed paths in its V9-graph resemble the computation of a sorting 
algorithm called Bubblesort. Its critical simplices are integer sequences with all 
entries positive. 

The second vector field is defined on K (N, 1), and it has only two critical 
simplices [] and [1], the same as the Eilenberg-MacLane field Veml- We call it 
the bit- chipping field and denote it by l^ch- 

Let us remark that one can consider composition of vector fields in a more 
general and more flexible setting, as is done in RSI 1 . but for our purposes, 
the simple notion above suffices. 



4.1 The bubblesort field 



Translating positive sequences to sorted sequences. To define the vec- 
tor field Vb s , it is convenient to consider a different representation of the sim- 
plices of K(Z, 1). Namely, we represent a fc-dimensional simplex a = [a%\ ■ ■ ■ \a^\ 
by a (k + l)-tuple (bo, b\, . . . , bk), where bo G % can be chosen arbitrarily and 
bi := bi-i + cij, i = 1,2, ... ,k. Thus, each a is represented as an equivalence 
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class of (k + l)-tuples of integers, where two (k + l)-tuples are equivalent if 
their difference is of the form (a, a, . . . , a) (all components equal). We denote 
the equivalence class of (bo, . . . ,bk) by [bo, ■ ■ ■ ,bk\. 

This correspondence between simplices of the form [a± \ ■ ■ ■ |etfc] and equiva- 
lence classes of (k + l)-tuples is obviously bijective. Nondegenerate simplices 
[ai|---|afc], i.e., those with no zero component, translate to [&o, •■■,&&] with 
bi-i ^bi, i = l,2,..., k. 

A (nondegenerate) simplex from K (N, 1) corresponds to [bo,...,bk] with 
strictly increasing components, i.e., 60 < &i < ■ ■ ■ < H- The face operators 
become extremely simple in this notation: di corresponds to deleting the ith 
component. 

The field. As was already announced, the critical simplices of Vb s are the 
[&o, • • • , bk] with bo < • • • < bk- If o = [bo, ■ ■ ■ , bk] is not critical, we look at 
the smallest i such that bg > bi + i; let us call it the leading index of a. Let 
us write v = be and u = 6^+1. We consider the maximal contiguous segment 
in the sequence 60,61,. . . starting at the £th position and containing only v's 
and u's; formally, we take the largest m > i + 1 such that bi € {u, v} for all 
i = I, t + 1, . . . , m, and either 6 m+ i {u, v } or m = k. We call b£, 6^+i, . . . ,b m 
the leading alternating segment of a (indeed, there can be no two consecutive 
u's or v's, since this would mean that a is degenerate), and we denote it by 
LAS(cr). 

Then we let a be a source if LAS(cr) ends with u, i.e., b m = u, and otherwise, 
a is a target. For a source a, still with u, v, m as above, we set 

T = Vb s (o-) := [6 , . . .,b m ,v,b m+1 , . . .,b k ], (2) 

i.e., V^s inserts another v just after LAS (a). 

With r = Vb s (cr) as in the just given definition, we have a = d m+ iT, and 
m + 1 is easily seen to be the only index i with a = d{T (thus, a is a regular 
face of t). Moreover, a can be uniquely reconstructed from r (delete the last 
element of LAS(r)), and so Vb s is indeed a discrete vector field. 

Lemma 4.1. Let a, r, I, m, u, v be as above. The only i 's for which a' := diT is 
a source simplex distinct from a are: 

(A) i = I, provided that m = 1+ 1 (i.e., LAS(cr) has length 2); then a 1 is 
obtained from a by swapping bi and bi + \ . 

(B) i = m + 1, provided that m + 2 < k, b m+ 2 = u, and LAS(o"') ends with u, 
where a 1 = d m +ir = [bo, . . .,b m ,v,b m+ 2 = u,b m +3, •••,&&]■ 

Proof. This is a straightforward case analysis. If a' should be a source simplex, 
then di must change LAS(r). It cannot delete elements from the middle of 
LAS(r), since the result would be degenerate, and it cannot delete the final v 
(since this was inserted by Vbs)- Thus, the only possibilities are de, which deletes 
the beginning of LAS(r), and d m+ \, which may make the leading alternating 
segment longer, if deleting b m+ \ connects LAS(r) to another u. 

The second possibility, <9 m +i, leads to (B). For the first possibility, d^, it 
remains to observe that if LAS(cr) starts with v,u,v,u, then LAS(c^r) starts 
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with v,u and ends with v (the one inserted by Vb s ), and so c^r is a target 
simplex. Thus, the situation must be as in (A). □ 

Corollary 4.2. The vector field Vt>s is admissible and polynomially bounded. 

Proof. First we observe that if we follow some directed path in the V9-graph, 



then once we make a step as in (B) of Lemma 4.1 we can never make a step 
of type (A) anymore, since the step of type (B) leaves us with the leading 
alternating segment of length greater than 2, where steps of type (A) cannot 
be applied. It is also clear that one can only make less than k consecutive steps 
of type (B), and then gets stuck. Thus, steps of type (B) can be ignored. 

A step of type (A), which swaps two consecutive elements from descending 
order to ascending, decreases the number of inversions in the current a, i.e., 
pairs with i < j and > bj. Thus, there is no directed cycle and Vt, s is 
admissible. 

Polynomial boundedness is now obvious: all simplices reachable from a = 
[bo, ... , bk] are made of integers in the set {bo, . . . , bp.}, and thus their number 
is bounded by a function of k (i.e., a constant from our point of view). A more 
careful analysis shows that at most 0(k 2 ) simplices are reachable from any k- 
dimensional a, and thus the bound for the bubblesort field is polynomial even 
as a function of k. □ 

4.2 The bit-chipping field 

Here we return to the "bar" notation [ai|a2| • • • |ofe], and we will consider only 
simplices of K(N, 1), which means a% > 1 for all i. 

The anatomy of a simplex. Let a = [ai|«2| • • • l a fc] he a nondegenerate 
fc-simplex of K(N, 1). We introduce the following terminology. 

• Let p = p(o) G {0, 1, . . . , k} be the largest index such that a\, . . . , a p are 
all powers of 2 and a\ < a2 < • • • < a p . The sequence ai|a2| • • • \a p is 
called the nondecreasing dyadic part of a. If 1 < p < k and a p > a p+ ±, 
then p is called the peak of a; otherwise, a has no peak. 

• Let q = q(a) G {0, 1, . . . , k} be the largest index such that a±, . . . , a q are 
all powers of 2 (thus, q > p). The sequence ai|ct2|-- - \a q is called the 
dyadic part of a. If q = k, then a is called fully dyadic. If, on the other 
hand, q < k, then q + 1 is the breakpoint of a and a g +i is the breakpoint 
value of a (which is not a power of 2). The sequence a q+ 2\a q+ ^ \ ■ ■ ■ \ak is 
the right part of a. 



s peak > breakpoint s breakpoint 

y v y 



Here are two concrete examples: 

t 

[ 2 2 | 2 3 | 2 6 | 2 6 | 2 2 | 2 5 | 7 | 4 | 3 | 12 ] [ 2 2 | 2 3 | 2 3 | 2 5 | 2 5 ^ | 38 | 5 | 8 | 11 ] 

nondecreasing right part nondecreasing dyadic right part 

dyadic part part = dyadic part 

^ ^ / 

dyadic part 
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The vector field. We define a vector field Vbch on K(N, 1). There are two 
types of source simplices. 

(a) The first type of source simplices are the simplices that are not fully 
dyadic and have no peak. Thus, all of the dyadic part is nondecreasing 
(i.e., p = q; we also admit p = q = 0) and the breakpoint value is larger 
than the last element of the dyadic part. Explicitly, they are of the form 

a = [2 41 \2 i2 \ ■■■ \ 2 i "\b\ a q+2 \ ■ ■ • | a k ] , 

2 n < 2*2 < ... < 2*9 < b. In this case we set 

V hch (a) =r:= [2 Jl | 2 i2 | • • • | 2 { " | lpow(6) | ltrim(6) | a q+2 \ ■ ■ • \ a k ) , (3) 

where lpow(&) is the largest power of 2 not exceeding b, and Itrim (6) := 
b — lpow(6). That is, r is obtained by splitting the breakpoint value b into 
two components, lpow(6) and ltrim(6); informally, we can think of this as 
"chipping off" the leading bit of b. 

We observe that each target simplex r as defined above has a peak, 
namely, p{r) = q(cr) + 1, and in particular, r has a nonempty dyadic 
part (but it may happen that the dyadic part of r is longer than the 
nondecreasing dyadic part, since Itrim (6) may be a power of two). 

(b) The second type of source simplices are the fully dyadic simplices a = 
[2 h | 2 i2 | • • • | 2 ik ] with 2 h <2 i2 <■■■< 2 1 "- 1 < 2 ik with i k > 1 (this last 
condition is important only for k = 1). In this case we set 

T = V bck ((r) := [2 h | 2' 12 | • • • | 2 1 '*- 1 | 2 ifc - x | 2 ik ~ 1 ] ; (4) 

i.e., we split the last component of a into two equal halves. 

Lemma 4.3. This definition indeed yields a vector field, and the only critical 
simplices are [} and [1]. 

Proof. Let us consider an arbitrary simplex r. If it is not fully dyadic and 
is not a source simplex, then it has a peak, and thus it has the form r = 
[2 h | • • • | 2*" | Cp+i | • • • | cfc+i] with 2 h < ■ ■ ■ < 2*? > c p+1 . This equals V hch (a) 
for a = [2 n | • • • | 2 1 p~ 1 \ 2 l ? + c p+ \ \ c p+ 2 | • • • | Ck+i]- Thus, r is a target simplex 
and there is exactly one edge (a, r) 6 Vb c h- Moreover, we have a = d p r, while 
djT ^ a for j ^ p, so a is a regular face of r as needed. 

Next, if r is fully dyadic and has a peak p, i.e., r = [2 n | • • • | 2 lk + 1 ], 2 n < 
■ ■ ■ < 2 lp > 2 lp+1 , then r is again a target simplex with r = Vbch(c) for a = 
[2*i | ••• | 2*'f- 1 | 2*p + 2*p+! | 2Vfi | • • • | 2* fc +i] (here 2*f + 2 1 'j>+i is the breakpoint 
value). Again, j = p is the only index with djT = a. 

The last remaining case is a fully dyadic r with no peak, which must be 
nondecreasing. If it is not a source simplex, then either we have one of the cases 
[], [1], or k > 2 and the last two components of r are equal, which means that r 
is of the form Q and a can again be uniquely reconstructed from it. We have 
a = djT for the unique index j = d — 1. □ 
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Preparations for analyzing Vb c h- It will be convenient to work mainly with 
the target simplices. Thus, given a target simplex r, we let t-reach(r) C reach(r) 
be the set of all target simplices reachable from r. 

First we will classify all possible target simplices r' reachable from a given 
target simplex r by two steps in the V9-graph; in other words, the r' of the 
form Vbc^djr) for some j. This is a straightforward, if somewhat lengthy, case 
analysis. The subsequent proofs of admissibility and polynomial boundedness 
will use this classification. From the point of view of these proofs, one could 
probably make some shortcuts and merge some of the cases in the classifica- 
tion; however, we prefer to give a complete picture, in order to avoid possible 
omissions. Indeed, for several other candidate vector fields we tried, "most" 
cases apparently worked fine, but those fields failed in what seemed like minor 
details. 

Lemma 4.4. Let r = [ai|a2| • • • \a k ] be a k- dimensional target simplex. 
If t is not fully dyadic, we can write it in the form 

j- 2 ii | 2*2 | ... | 2 %> | 2^+1 | ... | 2*« | b | a q+2 | • ■ ■ | a k ] , 

where b is not a power of 2, 2 n < • • • < 2 1 p , p > 1, p < q < k — 1, and either 
2 1 p > 2 l f +1 (if p < q) or 2 % f > b (for p = q). Let t' be a target simplex of 
the form Vh c h(djT) for some j, where a = djT is a (k — 1)- dimensional source 
simplex. Then t' has one of the following forms: 

(A) If p = 1 and 2 %2 < ■ ■ ■ < 2*« < b, then we can have 

t' = [2 i2 | ■ ■ ■ | 2*« | Ipow(ft) | ltrim(6) | a q+2 \ ■ ■ ■ \ a k ] 

(we drop the first component and split b). Example: r = [2 2 |1|2|7], r' = 
[1|2|2 2 |3]. 

(B) If ij < for some j, 1 < j < p — 1, then we can have 

t' = [2 h | • • • | 2^- J | 2^ +1 | 2^ | 2^+ 2 | ■ ■ • | 2*« | b \ a q+2 \ ■ ■ ■ \ a k ] 

(the entries 2 % i and 2^ +1 are swapped). Example: r = [1|2 2 |2|7], r' = 
[2 2 |1|2|7]. 

(C) If q > p + 2, i p — 1 = i p+ i = i p+ 2 < i p +3 < ■ ■ ■ < i g , and 2 % i < b, then we 
can have 

t' = [2 h | ■ ■ ■ | 2** | 2** | 2^+ 3 | ■ ■ ■ | 2^| lpow(6) | ltrim(6) | a q+2 \ ■ ■ ■ \ a k ] 

(two components following the peak are merged and b is split). Example: 
t = [2|1|1|2|7], t' = [2|2|2|2 2 |3]. 

(D) If q > p + 2 and i p+2 > ip > ip+i, then we can have 

t' = [2 h | ■ • • | 2** | 2 i '>+ 2 | 2^ +1 | 2*p+ 3 | • ■ ■ | 2*« | b \ a q+2 \ ■ ■ ■ \ a k ] 

(the entries 2 t v+ 1 and 2 t p+ 2 are swapped). Example: r = [2|1|2 2 |7], r' = 
[2|2 2 |1|7]. 
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(E) If q = p + 1, b' = 2 tp+1 + b satisfies b' > 2 lp , and b' is not a power of 2, 
then we can have 

t' = [2 h | ■ ■ ■ | 2 ip | lpow(6') | ltrim(6') | a q+2 \ ■ ■ ■ \ a k ] . 

Example: r = [2 3 |2|7], r' = [2 3 |2 3 |1]. 

(F) // the situation is as in (E) except that b' = 2 % is a power of 2, then we 
can have 

t' = V hch {[2^\---\2 ip \2 i \a q+2 \---\a k \) 

(note that here we do not write out r' explicitly, since there are still several 
cases to distinguish depending on the right part of t, but we will not need 
to discuss them explicitly). Example: r = [2 3 |1|7|19], r' = [2 3 |2 3 |2 4 |3]. 

(G) If q = p < k — 2, V :=b + a q+ 2 > 2 % p , and b' is not a power of 2, then we 
can have 

t' = [2 h | ■ ■ ■ | 2** | lpow(6') | ltrim(6') | a q+3 \ ■ ■ ■ \ a k ] . 
Example: r = [2 3 |7|4], r' = [2 3 |2 3 |3]. 

(H) // the conditions are as in (G) except that b' = 2 l is a power of 2, then we 
can have 

r' = V hch {[2^\-- -{2^12' \a q+3 \---\a k ]) 

(as in (F), we need not write out t' explicitly). Example: r = [2 3 |7|1|7] ; 
t' = [2 3 |2 3 |2 2 |3]. 

(I) If q = p = k — 1 and either p = 1 or i v -\ < i p , then we can have 

T ' = |~2*i | • • • | 2 ip ~ 1 | 2* p ~ 1 | 2* p_1 ] 

Example: r = [2|2 3 |7], r' = [2|2 2 |2 2 ]. 

If t = [2 n | • • • \2 tk ] is fully dyadic, then either p < k (t has a peak), or 
p = k (t is nondecreasing) and i k -\ = i k - I n the peak case, we have the 
following possibilities for t' = V^^djT): 

(dA) If p = 1 and i 2 < i 3 < • • • < i k -\ < i k , we can have 

r = [2 i2 | • • • | 2 ik ~ 1 | 2 ik ~ l | 2* fc_1 ] 

(deleting the first entry of t and splitting the last). 

(dB) For 1 < j < p — 1 and ij < ij+i, r' can be obtained by swapping 2 l J and 

(dC) If ' i p — 1 = ip+\ = ip+2 < ip+3 < • • • < ik~i < ik, we can have 
t 1 = |^2* 1 | • • • | 2 ip | 2 ip | 2 ip+s | • • • | 2 ifc_1 | 2 ik ~^ | 2* fc— ^] 
(merging two equal entries and splitting the last). 
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(dD) For i p+ 2 > i p > i p +\, r' can be obtained from r by swapping 2' 1 p +1 and 

Finally, if a fully dyadic r has no peak, we have the possibility (dB) for r' and 
the following additional one: 

(dl) If k = 2 or ifc_ 2 < ik-i> then we can have 

j-i — |~2*i | • • • | 2* fc ~ 2 | 2 ifc_1_1 | 2 ifc ~ 1 ~ 1 ] 

(drop the last component and split the previous one). 

Proof. As was already mentioned, the proof is totally straightforward and could 
probably be left to the reader. Yet, since getting used to the definitions and 
notation probably needs some practice, we chose to present the proof. 

As in the lemma, we first consider r not fully dyadic. If a = djT is a source 
simplex, then it has no peak, and thus the operation dj has to "destroy" the 
peak of r in some way. In particular, we have j <p+l, for otherwise, the peak 
of r is also present in djT. We just need to discuss the values of j in this range. 

For j = 0, do removes the first coordinate, and this may destroy the peak 
only for p = 1. For p = 1, a is a source iff 2 1 p +1 < • • • < 2 % i < b (this condition 
is void for q = 1), and if this holds, then r' is as in (A). 

If 1 < j < p - 1, a = [2 h \ ■ ■ ■ |2^- 1 |2^' + 2^+ 1 |2^+ 2 | • • • |2*«|6| • • •]. In this 
case, if ij = ij+i, then 2 l i + 2 l '+ 1 is a power of two, a necessarily has a peak, 
and thus it is not a source. So ij < ij+i; then a is a source and 2*? + 2 l i +1 is 
the breakpoint value, and r' is as in (B). 

Next, we consider j = p. Here the pth component of a is 2 %v + 2 lp+1 (for 
q > p) or 2 l f + b (for p = q). In both of these cases the pth component is not a 
power of 2 (since p was the peak of r), hence p is the breakpoint of a, and so 
^bch(c) = t. Therefore, j = p does not contribute any r'. 

Finally, we need to discuss j = p + 1 . Here the sum of the two entries of r 
following the peak must greater or equal to 2 z p (and, in particular, p < k — 2), 
for otherwise, p would be a peak in a. We consider three cases, depending on 
how many of these two entries are powers of 2. 

First, if q > p + 2, then the peak is followed by 2*p +1 and 2' 1 p+ 2 in r. 
If 2*p+i + 2^+ 2 = 2\ then i p+1 = i p+2 = i p - 1. Then a begins with 
[2*i | . . . |2 J p|2*p|2 l J'+ 3 | • • • \2 l i\b\ ■ ■ ■ , and since it has no peak, the dyadic part 
is nondecreasing. Then r' is as in (C). If, on the other hand 2 1 p +1 +2*p+ 2 > 2\ 
then 2 1 p +1 + 2*p+ 2 is not a power of 2. Then r' is as in (D). 

Second, we can have q = p + 1 (still with j = p + 1). Then the entry of a 
following 2 1 p is b' = 2 1 p +1 + b, which has to be at least 2 1 p. If b' is not a power 
of two, then r' is as in (E), and otherwise, we get (F). 

Third, we can have q = p. If q < k — 2, then the pth entry of a is followed 
by b' := b + a q+ 2, which has to be at least 2 % p . If b' is not a power of two, then 
t' is as in (G), and otherwise, we get (H). 

There is still one remaining case for j = p+ 1, namely, when p = k — 1; then 
dj just deletes the last coordinate and a is fully dyadic. Then a is a source 
precisely when p = 1 or i p _i < i p , and we have r' as in (I). 
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It remains to consider the case of r = [2 n \ • • • \2 lk ] fully dyadic; thus, q = k. 
First we assume that t has a peak p < k — 1. Then most of the analysis as 
above applies. 

For j = 0, we get that 8qt is a source iff p = 1 and «2 < h < • • • < < 
and then we have r' as in (dA). 

For 1 < j < p— 1, arguing as in the not fully dyadic case above, for ij < ij + \ 
we get t' by swapping 2^ and 2 i J+ 1 as in (dB). The case j = p again brings 



no t'. 



For j = p + 1, we have essentially the first of the three cases of the analogous 
analysis for the not fully dyadic case (q = k > p + 2). For i p — 1 = i p+ \ = 
ip+2 < V+3 < • • • < z/c-i < we obtain (dC), and for i p+ 2 > i p > i p +i we get 
(dD) (a swap). 

Finally, we may have r without a peak, which means that r = 
[2*i | ...| 2 ^-2|2^-i| 2 ifc-i], i x < ... < i fc _! (see case (b) of the definition of V hch ). 
Here 9 and d k ^ bring no r' (since V hch (d T) = Vb ch (9 fc _ir) = r). For 
1 < < k — 2 and ij < ij+i, we get a r' by swapping 2* J and 2 lj+1 as in 
(dB). For j = k, dk drops the last component, and if i k ~2 < ik-i, we get a r' 
by splitting the last component as in (dl). 

□ 



Acyclicity. Given Lemma 4.4, admissibility of Vbch can be proved quickly. 
Here we will check only acyclicity of the V9-graph, since the non-existence of 
infinite paths will be a side-product of the proof of polynomial boundedness 
below. 

Lemma 4.5. The Vd-graph contains no directed cycle. 



Proof. If r' = V\> c h{djT) is obtained from r as in Lemma 4.4, then for r not 
fully dyadic, one of the following can happen: 

1. g(r') > q(r), i.e., the length of the dyadic part increases. This is always 
the case in (F), (G), (H), and (I), and it may also happen in (A) and (C). 

2. q{r') = q{r) and the breakpoint value decreases. This happens in (A) and 
(C) (unless q drops) and also in (E). The latter is not entirely obvious, 
since we need to check that ltrim(2 lp+1 + b) < b, but this holds since 
ltrim(2%>+ 1 + 6) < 2 i f+ 1 + b-2 i ", and 2^ > 2 i ?+ 1 . 

3. q(r') = q(r), the breakpoint value stays the same, and the dyadic part 
becomes lexicographically larger. This happens in (B) and (D), since the 
swaps move a larger component forward. 

If r is fully dyadic, then so is r', and either the sum of components of r' 
is smaller than that of r (cases (dA) and (dl)), or the sums of components are 
equal and r' is lexicographically larger than r (cases (dB), (dC), and (dD)). 

This implies that there can be no directed cycle. □ 



We remark that an alternative proof of Lemma 4.5 can go along the following 



lines: If r = [ai \ ■ ■ ■ \ak] is not fully dyadic, then it can be shown that either 
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ones(r') < ones(r), where ones(r) is the total number of l's in ai, . . . , written 
in binary, or ones(r') = ones(r) and the sequence (ix,..., i p ) is lexicographically 
(strictly) larger than (i' l7 . . . ,i' pl ), where 2 H \ ■ ■ ■ \2 lp is the dyadic nondecreasing 

part of t, and similarly for 2*i| • • • |2V and r'. 

Polynomial boundedness. Condition (PBV1), polynomial computability 
of the vector field, is clearly satisfied for Vb c h> and so we need to check (PVB2); 
i.e., we need a polynomial bound on the total encoding size of all simplices 
reachable from a given simplex a. Obviously, we can focus only on target 
simplices: it suffices to provide, for every target simplex f, a polynomial bound 

On Eret-reach(f) S ' Ze ( T ) in termS ° f size(f). 

Moreover, it is easy to see that neither the application of T4> c h nor the face 
operators c\ can increase the sum of the components of the simplex. Thus, 
size(r) < size(f) for every r G t-reach(f), and it is enough to bound the number 
of simplices in t-reach(f). 

Thus, let us fix a target simplex f and set n := size(f). Our goal is a 
polynomial bound, in terms of n, on 1 1- reach (f)|. 

First we observe that fully dyadic simplices are easily accounted for. In- 
deed, a fully dyadic simplex [2 n | • • • \ 2 lk ] G t-reach(f) is specified by i±, . . . , if. G 
{0,1,... ,n — 1}, and so there are at most n k such simplices. 

So we consider only the r G t-reach(f) that are not fully dyadic. Let us 
write f = [5i| • • ■ |etfc] and r = [2 ll | • • • |2 l « |6|a g+ 2| ■ • ■ Qfc], where q = q(r) is the 
length of the dyadic part and b is the breakpoint value. 

We would like to show that with f fixed, there are only polynomially many 
possibilities for r. First, as was noted above, the number of choices for the 
dyadic part of r is polynomially bounded. 

Second, it turns out that all of the right part of r is inherited from f, i.e., 
ai = ai for all i > q + 2. This "stability of the right part" is not hard to prove 



inductively using Lemma 4.4 , and it will be the first part of the key lemma 
below. 

Thus, the last thing to do is showing that there are only polynomially many 
possibilities for the breakpoint value b of r, and this is the most tricky part of 
the proof. We will distinguish two cases: if b = a q+ \, i.e., b is "inherited" from 
f, then we call r a raw simplex, and otherwise, r is processed. 

The following lemma shows that if r is processed, then its breakpoint value 
belongs to a certain inductively defined set, which is of polynomial size. In order 
that the proof goes through, we need to strengthen the inductive hypothesis: 
namely, we need that for a processed r, the breakpoint value is smaller than the 
maximum entry of the dyadic part. This will play a role only in a single case 



among those in Lemma 4.4 namely (E); while all the other cases are natural 



and straightforward, (E) seems to work only by a small miracle. 

Lemma 4.6 (Key lemma). Let r G t-reach(f) be as above. Then ai = en for all 
i > q + 2, i.e., the right part of r coincides with the corresponding segment of 
t . Moreover, if t is processed, then b < max(2 n , . . . , 2 lq ), and b G B q+ \, where 
the sets B\, . . . , are defined inductively as follows: 

• B\ = ltrim*(ai), where, for a positive integer a, we define ltrim*(a) = if 
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a is a power of 2, and ltrim*(a) = {Itrim(a)} U ltrim*(ltrim(a)) otherwise. 

• Bj + i = ltrim*^{aj + i, aj+aj + i}U{2*+aj + i : < i < n—l}L){b+dj + i : b G 
Bj}^j, where we extend ltrim*(.) to sets by ltrim*(A) := |J agA ltrim*(a). 
Proof. It suffices to prove that if r is as claimed in the lemma, then t' = 



Vbch(djr) as in Lemma 4.4 has this form as well (moreover, we may assume that 



4.4 



t' is not fully dyadic). We need to consider the cases (A)-(I) in Lemma 
but we can right away settle (I), where t' is fully dyadic, as well as (B) and 

(D) , which only permute the dyadic part. This leaves us with cases (A), (C), 

(E) , (F), (G), and (H). 

First let r be raw, with b = a q+ i. In cases (A) and (C) r' contains lpow(a,j+2) 
followed by b' := ltrim(aq+2)i & t the (q + l)st position. If b' is a power of 2, then 
t' is raw, and otherwise, we have b' G B q+ \ and b' < lpow(a 9+ i); the latter is 
the required entry larger than b' in the dyadic part. Hence r' is a processed 
simplex as claimed in the lemma. 

In (E) and (G), we have a situation similar to the one just discussed, except 
that b' = a q+ i + 2 l for some i < n in (E), and b' = d q+ \ + a q +2 in (G). Moreover, 
in (E), b' is at position q + 1, while in (G) it is at position q + 2. Again we find 
that t' is a processed simplex of the claimed form. In cases (F) and (H), we 
either get r' fully dyadic, or the breakpoint value of r' is ltrim(a g / +1 ) for some 
q 1 > q + 1, preceded by lpow(o g ' + i). Then r' is a processed simplex as in the 
lemma as well, and the discussion of a raw r is finished. 

Now let r be processed, with b G B q+ i, b < max(2 n , . . . , 2 lq ). In cases 
(A) and (C) r may be raw, which is fine, or processed with breakpoint value 
ltrim(6), which lies in B q+ i, since B q+ \ is closed under ltrim(.). 

Case (E) is, in a sense, the most sophisticated, and it is here where the 
inductive hypothesis b < max(2 n , . . . , 2 tq ) is crucial. In the setting of (E), 2 lp 
is the maximum of the dyadic part of r, and so 2 %v > b. Let b' = b + 2 lp+1 , 
where 2*p +1 < 2*f; by the conditions in case (E), we have b' > 2 % ? . 

We claim that ltrim(6 / ) G ltrim*(6) (this will show that b' G B q+ \ and thus r' 
is as required). To check this, let us write, for brevity, u = i p and v = ip+i, and 
let (3 u -if3 u -2 ■ ■ ■ A) be the binary notation for b, i.e., b = Ylt=o Pi^i Pi e {0> !}• 
Since 2 U - 2 V < b < 2 U , we have = • • • = f3 v = 1. Then b' in binary is 

1000 • • • 0/?t,_ i/3„_2 ■ • • A), and so Itrim (6') can be obtained from b by iterating 
ltrim(.). Thus, b' G B q+ \ indeed. 

The consideration in cases (F) and (H) is the same as the one for f raw. 

The last case to consider is (G). Here the dyadic part of t' is longer than 
that of t. By induction, we have b G B q+ i, and so ltrim(& + a q+ 2) G B q+ 2 by 
the definition of B q +2 (or it is a power of 2, in which case r 1 is raw). As in 
the previous case, the entry lpow(6 + a g+ 2) supplies the power of 2 greater than 
ltrim(6 + a q+ 2), as required for the induction. The lemma is proved. □ 

Corollary 4.7. For f as in Lemma 4-6, we have |t-reach(f)| = 0(n 2k ), with 



the implicit constant depending on k. 

Proof. For each r G t-reach(f), we have at most n k choices for the dyadic part 
(which includes fixing q, the length of the dyadic part). A raw r is already 
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determined by f and by the dyadic part, while for r processed, we also need to 
specify b. 

The definition of Bj gives \B\\ < n and < 3n + n 2 + which 

yields \Bj\ = 0(n J ), and the corollary follows. □ 



Remark. A more careful (and more complicated) analysis should probably 
give 0{n k ) instead of 0(n 2k ) in Corollary 4.7 However, as we will now indicate, 



our vector field is not much better; there can indeed be about n k reachable 
simplices in t- reach (f). 

To see this, let us take n that is an integer multiple of k 2 , i.e., n = k 2 £, and 
let us consider a source simplex a = [d\\ ■ ■ ■ \dk], where di := (2 £ — 1)2^ -1 ^, 
i = 1,2, ... ,k. Put differently, if we think of the binary encoding of each di 
as consisting of k blocks of £ bits each (thus, di has at most n/k bits and 
size(cr) < n), then di has l's in the ith block and O's elsewhere. It can be shown 
that each simplex a = [ai \ ■ ■ ■ \ak], where a« has exactly one 1 in the ith block 
and O's everywhere else, belongs to reach (a). Since for each i, the position of 
the single 1 in m can be chosen in t ways, we have | reach(a)| > £ k = (n/k 2 ) k . 

It would be interesting to see if one could reach a significantly better bound 
with a different vector field, or if there is perhaps a good lower bound valid for 
every vector field. 
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